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A generalized CP symmetry for leptons is presented where CP transformations are part of an Sa 
symmetry that connects different families. We study its implications for lepton mixings in a gauge 
model realization of the idea using type II seesaw for neutrino masses. The model predicts maximal 
atmospheric mixing, nonzero #13 and maximal Dirac phase Sd — ±f • 



I. INTRODUCTION 



^vq The recent measurement of the leptonic mixing angle 613 in experiments searching for the oscillations of electron 

t-H antincutrinos emitted from reactors |T] and from the accelerator based experiments with muon neutrino beam [2] 
has generated considerable excitement in the field of neutrino physics. Taken together with already measured solar 
(N angle #12 and atmospheric angle #23, this almost completes the CP-conserving part of the lepton mixing matrix, 
O , under the assumption that there are no sterile neutrinos. This narrows the focus of the field to three remaining 
unknowns of neutrino flavor physics: (i) Dirac versus Majorana nature of the neutrino masses, (ii) mass hierarchy 
among them — namely, normal versus inverted — and (iii) leptonic CP-violating phases. The last item has two parts 
to it: Dirac phase, which is analogous to the CKM phase in quark sector, and Majorana phases, which are exclusive 
to the neutrino sector for Majorana neutrinos. The former can be measured in oscillation experiments whereas the 
latter may play a role in neutrinoless double beta decay searches. All these phases may play a role in understanding 
the origin of matter. 

On the theoretical side, despite such a vast amount of information, the nature of BSM physics responsible for 
neutrino flavor properties remains largely unknown and is the subject of extensive investigation. There are two 
generic approaches: one based on symmetries in the lepton sector leaving the quarks aside and a second one based on 
grand unified theories where both quarks and leptons are considered together. 

The quark-lepton unified grand-unification-based approach not only provides a very natural embedding of the seesaw 
mechanism to explain small neutrino masses but also, in a very economical class of renormalizable SO (10) models, 
turns out to be very predictive. Indeed, the recently measured value of #13 agrees with predictions made for this 
parameter in a minimal model of this type in 2003-2005 [3]. While this agreement is impressive, until there is some 
other evidence directly connecting the grand unification properties to seesaw physics (e.g., B-L violation as in 0]), 
qq one cannot test the GUT seesaw approach. 

The symmetry approach, on the other hand, derives its appeal from the fact that two of the observed neutrino 
mixing angles #23 (atmospheric) and #12 (solar) are close to values that look like group theoretical numbers, and 
find easy explanation in terms of simple discrete family symmetry-based models. For example, the observed near- 
£N) maximal atmospheric mixing can be easily understood if, in a basis where charged lepton masses are diagonal (to 
t-H be called "flavor basis" from here on), the Majorana neutrino mass matrix satisfies the Z 2 \i-r symmetry [5]. The 
simple versions of this symmetry, however, predict vanishing #13, a result which is contradicted by recent reactor [T] 
and accelerator experiments |2j . There is vast literature on the corrections to \x-r symmetry that come either from 
allowing general forms for the charged lepton matrix or from changing the neutrino mass matrix itself or combining 
simple ji-T symmetry with simultaneous CP conjugation [HI [7]. All these cases lead to nonzero #13. Many such models 
are also now ruled out since they predict values of #13 much smaller than the measured value. If, in addition to 
maximal atmospheric mixing, we consider the value of the solar angle tan #12 ~ we obtain the so-called tri- 
bimaximal mixing [5] , and it suggests more complicated groups such as Z 2 x Z 2 9 or S3 [TU] or A4 [IT] , but some of 
them also imply that #13 is zero or small after charged lepton corrections are taken into account and are not anymore 
phenomenologically viable. Thus the measurement of #13 has had a great impact on neutrino model building. 
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The discovery of large #13, however, does not rule out the generic symmetry approach and many examples have 
been discussed where new symmetries do allow for a large nonzero #13 |12H15j . We discuss one such approach in this 
paper which not only has the virtue of allowing large 6>i 3 but also predicts all the leptonic CP phases. The approach 
is somewhat different from many papers in the sense that we use a generalized definition of CP transformation among 
leptons |16j embedded in an S4 lepton family symmetry. We will call this new symmetry S4 symmetry. We present a 
gauge model for leptons invariant under this symmetry which not only accommodates a large #13 but also predicts a 
maximal 6*23 and maximal Dirac CP phase, i.e., = ±5. The maximal #23 is still consistent with the latest global 
analysis [TTJ, [TH] , although there are indications that it may be smaller [T7] . 

This paper is organized as follows: In Sec.JllJ we present the §4 model and the generalized CP transformation used 
in it; in Sec. |III| we present the various predictions of the model. In Scc. |IV| we give some comments and conclude 
with a summary of the results. In an appendix, we discuss the representations of the S4 symmetry that we use in the 
paper. 



II. MODEL 



Our model is based on the standard model gauge group SU(2)l <X> U(l)y with the usual assignment for leptons. 
Namely, the left-handed leptons Li transform as SU(2)l doublets with Y = — 1 and the right-handed charged leptons 
h ( = Ur) transform as singlets with Y = —2. The charged leptons gain masses through the Yukawa interactions 
with three Higgs doublets fa ~ (2,1), i = 1,2,3. Neutrino masses and mixing are generated through type II seesaw 
mechanism |19j which requires the introduction of Y = 2 SU{2) l triplets. In order to implement the symmetry in our 
model, we introduce four SM triplets, Ao and Aj ~ (3,2), i = 1,2,3, whose neutral members acquire small vacuum 
expectation values (vevs), induced by trilinear couplings of the form cf)(f>A'. We assume only three families of leptons 
and no singlet sterile neutrinos. 

We assume the theory to be invariant under a flavor symmetry acting in the horizontal space of the replicated 
fields. The chosen group is isomorphic to S4, but will contain generalized CP transformations (GCP) defined below; 
we denote this group by S4. Note that the group S4 has been pointed out as the group for tri-bimaximal mixing 20J, 
although some subgroup of it may turn out to be just accidental [HI [H]. The action of £4 on complex fields will be 
nontrivial. It is constructed as a subgroup of S4 <8> (CP) as follows. We remind the reader that S4 has generators S 
and T which satisfy the properties S 4 = T 3 = 1 and ST 2 S = T. 

Let us consider the (faithful) three-dimensional representation 3 of S4 generated by [55] 
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3:5= 0-1, T= 1 . (1) 
\0 1 0/ \o 1 0/ 

For complex fields, we can adjoin the usual CP transformation, denoted by the operator CP, to obtain S4 ® (CP). 
Note that S4 transformations and the CP transformation commute because all representations of S4 are real. We 
then extract the subgroup of S4 (g> (CP) generated by 



CP, T= 1 . (2) 





Notice the charge conjugation part in S is trivial for real fields. This group is isomorphic to S4 after we factor 
the subgroup generated by CP 2 = —1 for fermions. Such a factor group is §4. We keep the notation 3 for the 
representation generated by ^ . The other representations of S4 should be constructed in a similar manner from the 
representations 3', 2, l', 1 of S4. It is important to point out that S is a nontrivial GCP transformation that does not 
reduce to the usual CP transformation by basis change [TB] . 

Let us list the irreducible representations (irreps) of S4, constructed from the irreps 1,1', 2, 3, 3' of S4. They are 
led to peculiar representations of S4 when complex fields are considered: the real irreps 1 and l' (3 and 3') are 
interwoven in one equivalent (complex) representation 1 (3) whereas 2 splits into two inequivalent complex one- 
dimensional representations which we denote by l w and l u 2 ; see Appendix [X] for an explanation. They are quite 
similar to the representations 1', 1" of A4. 

We assign the representations of S4 as follows: 

Li — 3, Zi ~ 1, l 2 ~ l w , k ~ l^ 2 , <t>i ~ 3 ■ (3) 
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The fields assigned to the triplet representation Q transform explicitly as 



■ I . ;L .-«.<• ! : 

(4) 



Li(x) — > SijCL*(x) , Li(x) > TijLj(x) ; 
<f>i(x) -A Sij<j£(x) , <t>i(x) Tij4>j(x) , 



where x = (xq, — x) for x = (a?o,x) arises because of space inversion and C is the charge conjugation matrix. On the 
other hand, the right-handed lepton fields transform as 

h(x) — > Cl\{x) , h(x) — > h(x) ; 

h{x) -^Cl* 2 (x), ! 2 (i)^)wl 2 (i); (5) 
l 3 (x) — y Cl%{x) , i 3 (x) w 2 Z 3 (a;) . 
The Yukawa interactions for charged leptons invariant under these transformations is given by 



-£ y = yi (Zi0i + ^202 + L 3 4> 3 )h + y 2 {Li<t>i + ui"L 2 <j} 2 + ujL 3 4> 3 )1 2 
+ y 3 {L\(j>i + u>L 2 (/)2 + u 2 L 3 cj> 3 )l 3 + h.c, 

with the important restriction that all couplings yi are real due to invariance by S. 
When the neutral parts of the Higgs doublets acquire the vevs 



(6) 



^(M,l), (7) 



Mi = —= I 1 uj 2 w I diag(m e ,TO M ,m T ) . (8) 



the Lagrangian ^ gives rise to the charged lepton mass matrix 

'111 

1 LU 2 UJ 
: > \l CO u 2 

The correspondence is (m e ,m^,m T ) = v(y\,y 2 ,y 3 ) and we identify U* in ([8]) by defining 

1 I 1 1 l \ 

U u = -= 1 u uj 2 . (9) 
V 3 \1 w / 

We can see that MiMj has circulant form [5] and it is invariant by T and any transposition of family indices composed 
with complex conjugation (CP transformation), i.e., an 63 subgroup of S4. The matrix ^ is identical to the one 
obtained in A4 models. The potential for tpi is in fact the same as the general A4 invariant potential [TT], implying 
that A4 invariance leads automatically to S4 invariance for the potential of three Higgs doublets. For that potential, 
it has been shown that ([7]) is a possible minimum [23] . 

To generate neutrino masses, we introduce four Higgs triplets transforming under 54 as 

A ~l, A, ~3. (10) 

The S*4-invariant Lagrangian is then 

- C = i/oIfeAoL, + fx (IfeAiLa + ZfeA 2 Li + LfeA 3 L 2 ) + h.c, (11) 

where /o, fi are also real due to S. 

Given the large vev hierarchy, we can assume the potential allows arbitrary vevs for the neutral components of 
A 0> Ai, 

(A( 0) )= Mo , (Af)= Wi . (12) 



The Lagrangian (11) then induces the neutrino mass matrix 

fa ] 

M u =\f a d \ , (13) 
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where a = /qUq, d = fiu\,e = /1U2, / = /1M3. Notice the tri-bimaximal limit corresponds to e = / = [8 . For real 
a, d and complex e = /*, the symmetry corresponding to 23-transposition and complex conjugation (corresponding 
to an element of §4) would remain unbroken in the theory as symmetries of M v and M;M7. This would lead to CP 
invariance and nonzero $13. In contrast, if e = / we would obtain #13 = 0. In our case, CP violation and #13 7^ are 
allowed because there is no relation between e and /. 



If we assume the vevs ( 12 ) are real, the neutrino mass matrix, in the basis where the charged-lepton mass matrix 
is diagonal, is given by 




(14) 



where x, y are real while z is in general complex; they are independent combinations of the four parameters a, d, e, / in 
(13). This matrix has the same form as in [24] . invariant by [it exchange composed with complex conjugation (called 
/^T-reflection in [H]), with additional constraints so that it depends only on four real parameters. It has been shown 
that this form of the mass matrix leads to maximal #23 and maximal CP violation [7], with #13 7^ 0. 

The lepton mixing matrix Vmns wu 1 be the matrix that diagonalizes ( 14 ) . It is experimentally known that Vmns is 
close to the tri-bimaximal mixing matrix 



C/tb 



1 1 1 

" V^6 vj3 ;j/2 
\ 75 73 73 / 



(15) 



Therefore we parametrize 



where U e is the matrix that diagonalizes 



Vmns = U TB diag(l, 1, i)U e 



M' = U J M V U = 




(16) 



(17) 



e+f 
V2 



and 



e-/ 

V2 



with a,d,e,f being the original real parameters in (13). This 



for U = C/*[/ TB diag(l,l,i), b 

mass matrix has the same form as in the A4 model of [H] but our definition differs from [15] in that ( 16 1 includes an 
additional factor of i in diag(l, Therefore, our case corresponds to taking c purely imaginary in 12J. However, 
this case was not considered there because it was focused on nonmaximal # 23 and both real and imaginary parts were 



allowed to vary. In contrast, real a, d, e, / in the matrix ( 13 ) and, consequently, maximal #23, are natural consequences 
of our choice of symmetry. 

We can assume c > and consider the case c < by replacing i with — i in (16). This means that the sign of the 
Dirac phase 5d = ±f is not predicted in this model. Note that c controls #13 7^0 (and CP violation) and therefore 
it must be nonzero. 

The limit b, c — > leads to the tri-bimaximal form as U e = 1. As c 7^ to guarantee #13 7^ 0, U e should deviate 
from the identity. That means M' v must be nearly diagonal, i.e., |c| <C |a|, Having four parameters to describe 
9 quantities, we have 5 predictions, some of which are independent of the values of a, b, c, d. This is a consequence 
of the specific form of the mass matrix (14), i.e., maximal #23 and maximal CP violation 7J. In our specific model, 



the Majorana phases are also fixed: one is maximal and the other is zero. Only normal mass hierarchy for neutrinos 
is allowed. The remaining 5 physical quantities — two angles # 12 , # 13 and three neutrino masses mi, 7712,7713 — are 
correlated as they depend only on four parameters a, 6, c, d as discussed in the next section. 

A few comments are in order before we proceed to present the detailed numerical analysis of the model. 

• It is worth noting that in our model the lightest two neutrino eigenstates are almost degenerate in mass and are 
about a factor of three lighter than the third eigenstate unlike most normal hierarchy models where mil mi ~ 0.2 



• The Higgs potential for doublet fields in our model is the same as in the A4 models discussed in [33] and it is easy 
to see from there that there is a range of parameters in the scalar self-couplings where the vacuum alignment of 
the doublet fields in our model is justified. 



5 



III. PREDICTIONS 



In the limit b, c — > 0, the neutrino masses, i.e., the absolute values of the eigenvalues of (17 1, are given by 

|a + <i|, m 2 = |a| , m 3 = \a — d\ . 



mi 



(18) 



We can choose a > 0. From Am 2 2 = 



possibility. The experimental information Am| 3 = m 
(p| as 



m\ > 0, we can see that d < 0, hence normal hierarchy is the only 



3> Am 2 2 allows us to eliminate the modulus symbols in 



mi = \d\ — a ■ 



m 2 



m 3 = a 



\d\ 



We then arrive at the sum rule 



m 3 — 2to 2 — mi = , 



(19) 
(20) 



which commonly arises in models with discrete flavor symmetries |25| . The difference here is that the sum rule ( 20 ) 
applies to the neutrino masses themselves without additional Majorana phases or signs. 

When we allow b, c ^ 0, the sum rule (20) is still exactly satisfied provided that b = ±c. This can be seen from the 
eigenvalues of 



8M' V = M' v - ol 3 , 

which has characteristic equation 

- det(SM' v - Al 3 ) = A 3 — (d 2 + b 2 + c 2 )A - d{b 2 - c 2 ) = . 



(21) 



(22) 



The eigenvalues of M' v can be obtained from the roots of (22) by adding a. 



For general b and c the sum rule ( 20 ) is only valid approximately. The violation of the sum rule is quantified by 

(23) 



b 
d 



c 
d 



which controls the deviation of the PMNS matrix ( 16 1 from the tri-bimaximal mixing ( 15 ). The characteristic equation 

is 

(24) 



(22 1 shows that neutrino masses depend, apart from a, only on two combinations of d, c, b which can be chosen as 



d' = \d\Jl^4^? c , 5 



£ 213/2 



We can see 5 quantifies the violation of the sum rule. 

We can seek approximate roots to (22) for \8\ <C 1, which leads to 

—mi = a — d'(l — |(5) , 
m 2 = a — d'5 , 
m 3 = a + d'(l + ^5) . 



(25) 



The result is valid up to terms of order S 2 (order e 4 ) multiplied by d' . These relations can be inverted to write a,d' ,S 
in terms of the masses. In particular, the deviation of the sum rule is given by 



m 3 - 2to 2 - m% = ~S (m 3 + m x ) . 



(26) 



The knowledge of Am| 3 and Amf 2 determines the parameters a,d' in terms of 5. In turn, 8 depends on eb and e c , 
which affect #i 2 and 6*i 3 . 

To see how the mixing angles #i 2 and #i 3 are affected by ej,, e c , we can perform an analysis similar to [12 , with the 
difference that we have real matrices in our case. The matrix U e quantifies the deviations of the lepton mixing matrix 
from the tri-bimaximal form. For M' ul given the eigenvalues (— rni,m 2 ,m 3 ) in (25), we can calculate the eigenvectors 
which make up U e . The first approximation leads to 



1 e b 
-£fe 1 e c 



(27) 







1 
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where the real parameters £&, e c were given in (23 1. Notice e c > for c > because d < 



To first order, $13 depends only on e c , while #12 depends on e^, as 



We can then approximate 



sin 013 



sin 012 



3s 



ie 2 
3 c ' 

1 

3 



2^2 



(28) 



13 



1\ 2 



'12 



(29) 



where s 2 3 = sin 2 6*13 



and 



,2 _ „;„2 
s 12 



sin 0i2 as usual. This is the amount of deviation for the sum rule (26). We can see 



that the data [T5] are compatible with w 0. 

Given the experimentally known values of Amj 2 , Am| 3 , 0i2, 023, #13, we can determine the values of the neutrino 

masses 



mi 



13.3 meV, m2~15.9meV, m3«52.1meV. 



We have used the best-fit values of Rcf. 18 . A more precise numerical study reveals that 



11.8 meV < mi < 13.6 meV. 

when the la range for the observables is allowed [15]; see figures below. 
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Analogously, we can see the deviation for the sum rule is small as j5 ^ 
the deviation as 



(30) 



(31) 



0.1. In fact, our numerical study quantifies 



7Tl3 — 2?7l2 — mi 

7713 + nil 



= 11% to 15%, 



(32) 



at the la interval. 

The remaining numerical study is summarized in two figures. In Fig.[l] we display the range of sin 2 0i3 against 
the lightest neutrino mass. In Fig. [2] we display the effective light neutrino contribution m ee to neutrinoless double 
beta decay. Even though the two light neutrinos are quite degenerate in mass and have masses near 12 meV, due 
to Majorana phase, the effective mass is at most 3 meV. For both graphics the points are generated numerically 
without the analytic approximations employed in the previous analyses. We only collect the points compatible with 
the observables within 1-a as shown in Ref. H51. 



0.030 



0.028 



0.026 



0.024 



0.022 




12.0 12.5 13.0 13.5 

mi (meV) 



FIG. 1: Variation of #13 as a function of the lightest neutrino mass. This scatter plot was generated using e c > 0. There are 
points with e c < as well, corresponding to flipping the sign of the Dirac phase. However they do not introduce any perceptible 
change. 
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12.0 12.5 13.0 13.5 

nil (mcV) 



FIG. 2: The effective neutrino mass measured in neutrinoless double beta decay as a function of the lightest neutrino mass. As 
in Fig.[l] we have chosen e c > here. 

IV. CONCLUSIONS 

Wc have presented a model for leptons based on generalized CP symmetries which transform one family to another, 
generating the non-abelian S4 symmetry when supplemented by some permutations of families. This flavor symmetry, 
denoted by 64, represents a new implementation of the S4 symmetry where generalized CP symmetries are part of 
the group. This implementation shares some common features with the widely used group A4. For example, S4 also 
possesses three inequivalent one-dimensional representations that is similar to A4 in model building. The presence 
of CP transformations, however, further restricts the parameters of the Lagrangian to be real. The restrictions 
imposed by the generalized CP transformations are such that, with the addition of another Higgs doublet, we could 
have easily built another variant of the model where left-handed and right-handed leptons are assigned to the same 
representation 3 of S4. This could help us to embed this type of model in more symmetric theories such as left-right 
models. Therefore, this class of symmetries containing generalized CP transformations presents interesting features 
which can be further explored for flavor model building. 

Our specific model predicts maximal atmospheric mixing angle and accommodates the observed #13 without any 
cancellation among the model parameters; it predicts normal hierarchy and maximal Dirac phase of ±90° in the 
leptonic sector and should be testable in near-future long-baseline neutrino oscillation experiments. An important 
feature of the model is that the two light neutrino mass eigenstates are nearly degenerate in mass. Although the 
individual light eigenstates are "heavy," i.e., near 12 meV or so, due to maximal Majorana phase their net contribution 
to neutrinoless double beta decay amplitude is very small. The model also predicts an approximate sum rule relation 
valid for the three neutrino masses, without any Majorana phase or sign. The validity of the approximate sum rule 
is around 12%. 
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Appendix A: Other representations of S4 

We show here how to obtain the representations l w and l u 2 of §4 from the irreducible re pre sentation (irrep) 2 of 
S4. The irreps of §4 are constructed from the irreps of S4 by the procedure explained in Sec. O for the representation 
3: extract the subgroup of S4 <g> (CP) generated by S and T in ([2|, instead of S,T, CP that generate S4 ® (CP). 

Let us first explain why 1 and l' of S4 generate the same representation 1 of S4. We know that, for S4, 1 is trivial 
but 1' changes sign by S. If we follow the recipe and construct the representation of S4 corresponding to 1 and 1' we 



would obtain 



1 : S = S- CP-> 1-CP, T->1, 
1': S = S- CP ->(-!)■ CP, T-S-l. 



(Al) 



We are using the generators Q of 3 of S4 as the group elements themselves, given that the representation is faithful. 
The CP transformation denoted by CP acts as usual. A fermion field ip(x) and a complex scalar field <f>(x) transform 
as 



ip(x) Cip*(x) 
<f>(x) 4>*{x) . 



(A2) 



Therefore, the representations (All are equivalent because if ip(x) transforms as 1 then ip'(x) = iip{x) transforms as 
l' (notice that the field must be complex) . This same reasoning leads to the equivalence of 3 and 3' when we go from 
S4 to S4. 

Let us see what happens to the representation 2 of S4, which is equivalent to the S4 — > S3 homomorphism. To 
preserve the structure of S4 ® (CP) for which CP commutes with the elements of S4, it is important to consider real 
representations of 2. We adopt a slightly different version of Ref. (32 Eq.(640), in second reference], 



D 2 (S) 



1 

-1 



D 2 (T) 



1 v(3 s 

2 2 
vd 1 

2 2/ 



(A3) 



Then S4® (CP) in this representation is generated by D 2 {S), D 2 (T) and CP acting as in ( A2 1. The subgroup S4 would 
be generated by 



D 2 (S) 



1 

-1 



CP, D 2 (T) 



1 v(3 s 

2 2 
V3 1 

2 2/ 



(A4) 



However, it is usually more convenient to work with the complex basis where T is diagonal. We change basis to 

'1 



where 



with the basis-change matrix 



D' 2 (S)= I) CP, D' 2 (T) = 
D' 2 {T) = Jft D 2 (T)X , 



w 
to 2 



X = 



A 1 



Now, D^S) in (A5I differs from 



1 

-1 



X ■ CP 



(A5) 



(A6) 



(A7) 



(A8) 



because X is complex and complex basis change acts differently for CP transformations. The correct transformation 
is 



D' 2 (S) = Xt 



1 

-1 



X* ■ CP. 



(A9) 



which leads to ( A5 ) . 



Equation ( |A5[ ) defines the representation of S4, derived from 2 of £4. Since both transformations which generate 
S4 do not mix the first and second components, they are essentially one dimensional (complex). They correspond 
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to the representations which we denoted by l u and l w a , corresponding to the action of ( A5 ) to the hrst and second 
components, respectively. Explicitly, for a fermion field ip(x) (chiral or not), we have 

l u :i/>(x)—i-Cip*(x), ip(x) wip(x) , 
l u i : iP(x) A Cty* (i) , tl>(x) A ^V(x) . 



U) 1 



If could ignore gauge quantum numbers, the representation \ u and l w 2 would be equivalent because if ^>(je) ~ 1 
then Cip*(x) ~ l w 2. Its real representation space is two dimensional. In particular, 1 and l' would correspond to 
CP-even and CP-odd combinations of fields which have no definite transformation properties under the gauge groups. 
It is important to emphasize that if we were considering the whole S4 ® (CP), the representation 2 would remain two 
dimensional (complex). 
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